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We describe an explicit construction of approximate Ginsparg- Wilson fermions for QCD. We use ingredients of 
perfect action origin, and further elements. The spectrum of the lattice Dirac operator reveals the quality of the 
approximation. We focus on f3 = 6 for optimisation. Such fermions are intended to be inserted into the overlap 
formula. Hence we also test the speed of convergence under polynomial evaluation of the overlap formula. 



1. Introduction 

If a lattice Dirac operator D obeys the 
Ginsparg- Wilson relation (GWR) 

{D x , y , 75 } = 2{D 75 RV) x , y , (1) 

then it reproduces correctly the physical proper- 
ties related to chirality In particular, a GW 
fcrmion has a chiral symmetry, which is lattice 
modified but exact ||. 

The kernel R must be local, and {^,75} 7^ 0. 
For simplicity we assume 

Rx.y = -^Sx, y (A* > 0), and D ] = J5D75. (2) 

Then the spectrum of D is situated on a circle in 
the complex plane: it is the circle through with 
centre and radius /i. 

Classically perfect fermion actions solve the 
GWR H and their scaling is excellent, but they 
are hard to construct. They were applied success- 
fully in the Schwinger model H,[| , and the tedious 
work in d = 4 is still in progress Q . 

There are, however, many more GW fermions: 
starting from almost any sensible lattice Dirac op- 
erator Do (local, free of doublers etc.) we can 
generate a solution to the GWR by means of the 
overlap formula 

D ov = Jl + -==), A = D -ii. (3) 

(The allowed range of fi > will be commented 
on below). If we start from the Wilson fermion, 
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Do = Dyy, we arrive at the Neuberger fermion J7| . 
However, here we are interested in improvements 
in various respects due to a better choice of Do 

§• 

Our guide-line for the search of a better Do is 
the following observation: if Do is a GW opera- 
tor already (for some fixed GW kernel R), then 
it coincides with the resulting overlap operator, 
D ov = Dq. In real life we are now going to 
construct an approximate GW operator within 
a short range, and we use it as Dq- Hence the 
notorious square root in eq. ([|) is approximately 
constant, V A^A ss /1, and therefore 

Dov « Do . (4) 

For the free fermion, the perfect action can 
be constructed and parametrised explicitly 
though it involves an infinite number of couplings. 
Still this fcrmion formulation is local, because the 
absolute value of these couplings decays exponen- 
tially with the distance between ip and ip. 

We construct Do from the couplings of a trun- 
cated perfect free fermion. The truncation is done 
so that only couplings inside a unit hypercube 
survive, hence we arrive at a hypercube fermion 
(HF). Then we gauge this HF (we distribute the 
coupling over the link paths) by using just very 
few parameters, in such a way that the violation 
of the GWR is minimised. The criterion for this 
violation is proximity of the fermion spectrum to 
the unit GW circle, which corresponds to fi = 1. 

Even if we only perform "minimal gauging" 
(shortest lattice paths only, see below), the re- 
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suiting operator Dhf is promising for scaling and 
it is almost rotationally invariant jl(J — due to 
the perfect action background — although it is 
short-ranged. 

If we now insert Do = Dhf in the overlap 
formula, we expect a number of virtues for D ov , 
which are all based on relation (|J): 

• D ov is promising for a good scaling be- 
haviour and approximate rotational invari- 
ance, properties which are inherited from 
D . 

• We expect D ov to have a high level of local- 
ity (a fast exponential decay of the correla- 
tions) , since Do is ultralocal and long-range 
couplings can be turned on just slightly by 
the overlap formula. 

• The iterative transition from Dq to D ov is 
fast, since we already start off in the right 
vicinity. 

This is in contrast to the Neuberger fermion, 
which does a rather poor job with respect to 
all these properties. The above predictions have 
been tested and confirmed in a comprehensive 
study of the 2-flavour Schwinger model We 
used a 2d HF, which had (amazingly) a similar 
scaling quality as the very mildly truncated, clas- 
sically perfect fermion of Ref . [Q . The superiority 
of the resulting D ov (improved overlap fermion) 
over the Neuberger fermion (or standard over- 
lap fermion) was striking in all the respects listed 
above. 

We are now carrying on this program to 
(quenched) QCD |13|. Locality properties of free 
4d improved overlap fermions have already been 
discussed in Ref. Q, and first results for QCD 
have been reported in Ref. Jl2| . The same concept 
has also been adopted in Ref. (m) , which presents 
some results with respect to the last property 
(speed of the transition), based on a very simple 
choice of Do- On the other hand, sophisticated 
and complicated approximate GW fermions have 
been constructed in Ref. fill . 

Following the same lines, one may also work on 
improved domain wall fermions by optimising the 
choice of the 4d operator that corresponds to D 
S. Recent proposals focus on accelerating the 



transition to exact chirality, which is expressed 
here in terms of the required extent of the extra 
dimension ]lq] . 

2. Construction of Dq 

We now describe explicitly the construction of 
our HF, which approximates a GW fermion, and 
which is also promising in other respects. We 
proceed gradually in a sequence of steps. 

2.1. Truncated perfect free fermion 

As we mentioned earlier, perfect actions for free 
fermions can be constructed analytically, where 
the term R occurs in the renormalisation group 
transformation term. The locality is optimal for 



Rx. 



\5 X . V (i.e. fi = 1) 



which is the stan- 



dard GW kernel, and which we are going to use 
in the following. 

We write the Dirac operator of the free perfect 
fermion as D(x - y) = p^x - y)^^ + X(x - y), 
where we denote p M , A as the vector term and 
the scalar term, respectively. We truncate both 
terms by imposing periodic boundary conditions, 
so that we arrive at a free HF. The explicit cou- 
plings inside a unit hypercube are given in Ref. 
p0|| , Table 1. The scaling of this free HF is excel- 
lent jl0| , and its spectrum is very close to a GW 
unit circle (81, hence it is a good approximation 
to a GW fermion. 

The resulting free overlap HF has the expected 
virtues, in particular good scaling and locality j^j. 

2.2. Minimal gauging 

We now proceed to quenched QCD, and we al- 
ways use the plaquette gauge action. Tests with 
improved gauge actions are on the way. 

Our point of departure is a minimal gauging 
of the HF: the free couplings are attached to the 
shortest lattice paths only, in equal parts where 
several shortest paths coexist. This is the sim- 
plest version of a HF in gauge theory, but the 
harmony of the truncated free couplings leads 
to a beautiful pion dispersion relation, even at 
strong coupling, see Ref. |l0| , Fig. 10. On the 
other hand, this kind of gauging implies a strong 



They seam to be indeed useful for chirality, see last ref- 
erence in 1151 . 
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mass renormalisation — comparable to the Wil- 
son fermion. The reason is that all paths are sup- 
pressed by the gauge field (compared to the free 
case), whereas the coupling between if and ip on 
the same site in unaltered. For instance, at (3 = 5 
we obtain a "pion mass" of 3.0 jnj. 

This effect is also visible from the fermionic 
spectrum: again at (3 = 5 the right arc still fol- 
lows closely the GW circle — unlike the Wilson 
fermion spectrum, which extends to large real 
parts. However, the physically crucial left arc is 
absent in both cases, see Fig. [l]. 



p=5: Wilson fermion 
hypercube fermion 




Figure 1. The fermionic spectrum for the mini- 
mally gauged HF and for the Wilson fermion, in 
a typical configuration at j3 = 5. 

In such a situation, the overlap formula - 
which is supposed to provide a projection of the 
spectrum onto a GW circle — is not applicable 
any more: there is no suitable value of [i, which is 
the centre of the circle to be mapped on. Wher- 
ever we choose it, we would be confronted fre- 
quently with too many mappings to the left (re- 
vival of the doubling problem) or to the right arc 
(mass renormalisation is back) or both. 

Now it is of interest where the applicability of 
the overlap formula sets in. At j3 m 5.4 the cou- 
pling is still too strong, whereas at (3 > 5.6 we are 
statistically on safe grounds for a simple operator 
Do 0. Q' 



3 Note that this is a situation where it could be misguiding 
to consider only the spectrum of A' A, as it is sometimes 
done in the literature. 

4 If we deform a configuration continuously into another 
topological sector, then at least one eigenvalue has to cross 
the centre of the circle, since we define the topological 



So far, we show the full fermionic spectra on 
small 4 4 lattices. From our experience, this is 
sufficient to obtain the qualitatively correct im- 
pression about the performance of a lattice Dirac 
operator. The fact that the arc close to zero tends 
to be absent is not due to the fermion formula- 
tion, but solely due to the small lattice. Later we 
will show this arc specifically also on 8 4 lattices, 
in order to provide a more complete picture. 

2.3. Criticality through link amplification 

We now go beyond minimal gauging by intro- 
ducing a few extra parameters and optimising 
them so that the GWR violation is minimised at 
[3 = 6. 

First we attach a amplification factor l/u 
(u < 1) to each link. The idea is to compensate 
the mean link suppression due to the gauge field. 
At (3 — 6 we obtain criticality at u ~ 0.8. This al- 
ready provides a decent approximation to a GW 
fermion, see Fig. 0. 
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Figure 2. HF spectra at j3 = 6 on a 4 lattice 
at critical link amplification l/u without fat links 
(a = 0) and with fat links (a = 0.3J. 

2.4. Fat links and vector term suppression 

We now introduce fat links by globally substi- 
tuting all links as 



link 



(1 — a) link + — stapl 



It turns out that for instance a = 0.3 is a useful 
value: it helps to move the upper (resp. lower) arc 



charge by the index theorem. What we achieve at /3 > 5.6 
is that this transition is very quick, i.e. for most config- 
urations occurring at that coupling strength, all (almost) 
real eigenvalues clearly belong to the left or the right arc. 
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closer to the circle, and it also pulls the eigenval- 
ues closer together, which are both desired effects, 
as illustrated in Fig. g. 

On the other hand, a sequence of experiments 
suggests that the clover term is not a very pow- 
erful tool to improve the GWR approximation. 
Hence we omit it, since we only want to include 
new ingredients if they really yield a significant 
progress. 

An exception from that rule are the factors that 
we attach to the links, since they require practi- 
cally no computational effort. Hence we can also 
introduce different link factors for the scalar term 
A and the vector term p M . The former is responsi- 
ble for the chiral limit, so we leave it unaltered at 
U = 0.8. However, we suppress the vector term a 
little, because it is responsible for the imaginary 
part of the fermionic spectra, which are a bit too 
large so far. So now we multiply in p^ only the 
links by v/u, v < 1. By tuning the value of v we 
can make the above spectra follow the shape of 
the GW circle. Still the fat link is useful to sup- 
press the radial fluctuations of the eigenvalues. 
We use again a = 0.3, and find the suitable vec- 
tor link suppression to be v = 0.92. This leads to 
the quite satisfactory spectrum shown in Fig. 
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a=0.3, u=0.76, v=0.92 
L=8 
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Figure 3. The HF on a 4 4 lattice with critical 
link amplification, fat links and a suppression of 
the vector term. We also show the "continuation" 
around on a 8 4 lattice with the same parameters, 
i.e. u = 0.76, a = 0.3, v = 0.92. 



3. Convergence to D ov 

In d — 4 the troublesome square root in the 
overlap formula (0) can only be approximated it- 
eratively. A simple approach is the expansion of 
the inverse square root in a series in A* A — p? . 
This converges for all eigenvalues only if we start 
from a good G W fermion approximation Dq , and 
then the convergence is fast, as we have shown in 
the Schwinger model 

We now turn our attention to the usual method 
which approximates the sign function e in 

D ov = n[l + 7se(H)] , H = Tj (D - p) = H f (5) 

by a polynomial in H. 

The histograms for typical spectra of Hw and 
Hhf at (i — 6 are shown in Fig. || (on top). 
For the HF (with a = 0.3, u = 0.76, v = 0.92) 
the spectrum is already nicely peaked near ±1, 
whereas the Wilson fermion yields a broad distri- 
bution. 




Figure 4. On top: The eigenvalue histograms for 
typical spectra of Hhf (grey) and Hw (black) at 
(3 = 6. Below: The same after re-scaling so that 
all | eigenvalues \ are < 1 . 

In order to make the polynomial approxima- 
tion applicable for all eigenvalues, we have to re- 
scale the spectra such that they are confined in 
the interval [—1,1], see Fig. || (below). This does 
not affect the spectrum of Hhf too much, and 
in particular there is a large gap around 0, where 
any polynomial approximation of the sign func- 
tion has its largest error. On the other hand, for 
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Figure 5. The spectra of approximate overlap 
fermions, where the sign function is replaced by a 
polynomial of degree 21. We use a configuration 
typical at (3 = 6 and start from Dq = Dyy ( stars ) 
resp. from Dq = Dhf (crosses). 

Hvv there is a significant eigenvalue density in the 
vicinity of 0, which shows again that here it is far 
more demanding to enforce convergence to a GW 
fermion. 

We illustrate the difference by using a linear 
combination of Chebyshev polynomials (as sug- 
gested in Ref. Jl7|] ) with maximal degree 21 to 
approximate the sign function in the overlap for- 
mula (|). For a typical configurations at (3 = 6 
(on a 4 lattice) this leads to the spectra shown in 
Fig. |5| for the Wilson fermion resp. the HF. We see 
that the latter is clearly superior. This gain might 
already compensate the computational overhead 
by a factor of O(10) in the application of Dhf- 

4. Conclusions 

A good approximation to a GW fermion in 
QCD at [3 = 6 has been constructed, using only 
10 independent parameters. This HF is ready 
to be inserted in the overlap formula (||). This 
works down to j3 « 5.6, and it leads to an ex- 
act GW fermion. There are good reasons to ex- 
pect that the resulting GW fermion is superior 
over the standard Neuberger fermion with respect 
to scaling, rotational symmetry and locality. We 



5 A restriction to degrees below fa 20 is not sensible, re- 
gardless what Do one uses, because the approximation to 
the sign function is too poor. This is also true for other 
polynomials that may be used for this purpose. 



tested the speed of convergence towards an over- 
lap fermion if the sign function in the overlap for- 
mula is replaced by Chebyshev polynomials. Also 
with this respect we found a significant gain if we 
start from Dhf instead of D\y- This effect will 
be discussed in detail in Ref. Jl3| . 
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